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Abstract 

We study the solutions of the nonstationary incompressible Navier-Stokes equations in 
M. d , d > 2 , of self-similar form u(x,t) — -fiU(-j=) , obtained from small and homogeneous 
initial data a(x) . We construct an explicit asymptotic formula relating the self-similar 
profile U (x) of the velocity field to its corresponding initial datum a{x) . 

1 Introduction 

In this paper we are concerned with the study of solutions of the elliptic problem 
i-lU-l( X -V)U-AU + (U-VU )+ VP = 

o: \v-c/ = o, 

where U = (Ui, . . . , Ud) is a vector field in M. d , d > 2, V = (d\, . . . ,da) , and P is a 
scalar function defined on R rf . Such system arises from the nonstationary Navier-Stokes 
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equations (NS), for an incompressible viscous fluid filling the whole Mr, when looking for a 
velocity field u(x,t) and pressure p(x,t) of forward self-similar form: u(x,t) = -^U(xj\fi) 



o 



and p(x,t) = jP(x/y/i) . An important motivation for studying the system (1) is that the 
corresponding self-similar velocity fields u(x, t) describe the asymptotic behavior at large 
^ ■ scales for a wide class of Navier-Stokes flows. Moreover, simple necessary and sufficient 

conditions for a solution of the Navier-Stokes equations to have an asymptotically self-similar 
profile for large t are available, see [16]. We refer to [4] and [13], for more explanations and 
further motivations. 

The problem that we address in the present paper is the study of the asymptotic behavior 
for \x\ — > oo for a large class of solutions to the system (1). 

The existence of nontrivial solutions of (1) has been known for more than sixty years. 
For example, in the three-dimensional case Landau observed that, putting an additional axi- 
symmetry condition one can construct, via ordinary differential equations methods, a one- 
parameter family (U, P) , smooth outside the origin, and satisfying (1) in the pointwise sense 
for x / (see, e.g., [1, p. 207]). 



'Keywords: asymptotic profiles, asymptotic behavior, far-field, large distance, selfsimilar, incompressible 
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Landau's solutions have the additional property that U is a homogeneous vector field of 
degree —1 and P is homogeneous of degree — 2 , in a such way that the corresponding solution 
(u,p) of (NS) turns out to be stationary. A uniqueness result by Sverak, [18] implies, on the 
other hand, that no other solution with these properties does exist in M 3 , other than Landau 
axi-symmetric ones. See also [12] for a detailed study of the asymptotic properties of these 
flows. 

The class of solutions to the system (1) is, however, much larger. Indeed, Giga and 
Miyakawa [10] proposed a general method, based on the analysis of the vorticity equation in 
Morrey spaces, for constructing nonstationary self-similar solutions of (NS). A more direct 
construction was later proposed by Cannone, Meyer, Planchon [5], [6], see also [13, Chapt. 
23]. Now we know that to obtain new solutions U of (1) we only have to choose vector fields 
a(x) in M. d , homogeneous of degree —1, and satisfying some mild smallness and regularity 
assumption on the sphere : the simplest example in M 3 is obtained taking a small e > 
and letting 

but a condition like a|§d-i € L°°(S rf_1 ) with small norm (or similar weaker conditions involv- 
ing the L d -norm or other Besov-type norms on the sphere) would be enough. The basic idea 
is that the Cauchy problem for Navier-Stokes can be solved, through the application of the 
contraction mapping theorem, in Banach spaces made of functions invariant under the natural 
scaling. The profile U of the self-similar solution u obtained in this way (i.e. U = u(x, 1) ) 
then solves the elliptic system (1). 

Regularity properties and unicity classes of those (small) self-similar solutions have been 
studied in different functional settings (see, e.g., [14], [9]) and are now quite well understood. 

On the other hand, probably because of the lack of known relations between the self-similar 
profile U and the datum a , even in the case of self-similar flows emanating from the simplest 
data, such as in (2), the problem of the asymptotic behavior of the solutions U obtained in 
this way has not been addressed, before, in the literature. 

The main purpose of this paper is to construct an explicit formula relating U(x) to a(x) , 
and valid asymptotically for |x| — > oo. 

We will also consider the more general problem of constructing asymptotic profiles as 
\x\ — > oo for (not necessarily self-similar) solutions u(x,t) of the Navier-Stokes equations with 
slow decay at infinity (tipically, \u(x,t)\ < C|x| _1 ). Our motivation for such generalization 
is that solutions with such type of decay have, in general, a non-self-similar asymptotic for 
large time. In fact, Cazenave, Dickstein and Weissler showed that their large time behavior 
can be much more chaotic than for the solutions described by Planchon [16]. As shown in 
[7], however, one can obtain some understanding on the large time behavior of these solutions 
from the analysis of their spatial behavior at infinity. 

2 Main results and methods 
2.1 Notations and functional spaces 

If Q = (Qj-h,k) an d B = (Bh^k) are, respectively, a three-order and a two-order tensor in the 
Euclidean space M. d , we denote by Q : B the vector field with components 

d 

(Q:B)j= Qj;h,kB h ,k, j = l,...,d. 

h,k=l 

Sometimes, in the proofs of our decay estimates, we will simply write QB instead of Q : B 
when all components of such vectors can be bounded by the same quantities. 
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We denote the gaussian function by 

gt ( x ) = (47rt)- d / 2 e- |x|2/(4i) , xeR d , t>0. 

As usual, we adopt the semi-group notation e tA a = g t * a for the solution of the heat system 
dtu = An, with u\t=o = a for an initial datum a defined on the whole M d . 

All the functions we deal with are supposed to be measurable. By definition, for any $ > 
and m G N, 

/ G £™ / G C m (M\{0}) and Ix^ 1 " 1 9 a / G L 00 (M a! ) Vq G N d , |a| < m. (3) 

We are especially interested in the case •& = 1. Indeed, the spaces -E™ contain homogeneous 
functions of degree —1 (and, in particular, the initial datum a(x) given by (2)). 

The non-homogeneous counterpart of Eff is the smaller space E™ , which is defined by 
the additional requirement that d a f G L°°(M. d ) for all \a\ < m. These spaces are equipped 
with their natural norm: 

ll/H = max sup \xf + \ a \\d a f(x)\, 
o \a\<m x< zM. d \{0} 

\\f\\ ET = max sup(l + Nr+H|^ /(x) |. 

\a\<m xeRd 

Our starting point is a classical result by Cannone, Meyer and Planchon about the con- 
struction of self-similar solutions of the Navier-Stokes equations 

dtu + V ■ (u <g) u) = Au — Vp 
< divu = x G R d , t > 0. 

u\ t =o = a, 

Even though their construction goes through under very general assumptions on the regularity 
of the initial data, here we are mainly interested in the following simple result: 

Theorem 1 (see [5], [6]) For all m G N there exist e, (3 > such that for all initial datum 
a G E™ homogeneous of degree —1, divergence-free and satisfying 

\\a\\^rn < e, (4) 

there exists a unique self-similar solution u(x,t) = -^u(^-^j of the Navier-Stokes system 

(written in the usual integral form, see (NS) below) starting from a, and such that \\U\\e™ < 
(3. Moreover, 

U(x) = e A a + C(|xr 2 ), as \x\ -> oo. (5) 

More precisely, Cannone, Meyer and Planchon prove that U(x) = e A a(x) + TZ(x) , where 
the remainder term satisfies 1Z G E^ 1 ■ Their result was stated in dimension three, but their 
proof easily adapts for all d > 2 . 

2.2 Main results 

Our main result shows that one can give a much more precise asymptotic formula between 
the asymptotic profile U (x) and the datum a(x) . It turns out that such asymptotic profile 
has a different structure in different space dimensions. 
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Theorem 2 Let a(x) he a homogeneous datum of degree — 1, such that a is smooth on the 
unit sphere S^ 1 and satisfying the smallness condition (4) for some m > 3. Let u(x,t) = 
■±=Jj(^-^J the self-similar solution constructed in Theorem 1. Then the following profiles 
hold: 

• If d = 2, we have as \x\ — > oo, 

U(x) = a(x) - log(|x|)^j^ + 0(|x|- 3 ), (6a) 

Here A = (A h>k ) is the 2x2 matrix given by A hjk = J sl (a h a k ) and Q(x) = (Q j; h,k( x )) > 
where the Qj-h,k are homogeneous polynomials of degree three (given by the explicit 
formula (12) below) 

• For d = 3, we have as \x\ — > oo, 

O(r) ■ B 

U(x) = a{x) + Aa{x) - PV ■ (a ® a) - j^" + 0(|xr 5 log |x|) , (6b) 

for a dx d constant real matrix B = (B^ k ) depending on a. Here P = Id — V(A) -1 div 
is the Leray-Hopf projector onto the divergence-free vector fields. 

• For d > 4, the far- field asymptotics reads, as \x\ — > oo, 

£/(x) = o(x) + Ao(s) - PV • (a O a) + 0(|x| -5 log |x|) . (6c) 



In Section 8 we will restate and prove this theorem in a more general form, removing the 
assumption that a is homogeneous. Such more general theorem will apply also for solutions 
u(x, t) of Navier-Stokes of non-self-similar form. On the other hand, we will not seek for the 
greatest generality about the regularity of the datum: even though there is a considerable 
interest in studying self-similar solutions emanating from rough data (see [13], [11]), in most 
of our statements we will assume that a G C^S* -1 ), which is of course non-optimal, but 
permits us to greatly simplify the presentation of our results and to better emphasize the 
main ideas. 

The method that we present in this paper would allow to compute, in principle, the 
asymptotics of U up to any order, when a is smooth on S d_1 . However, the higher order 
terms have quite complicated expressions. 

The functions Aa and PV ■ (a® a) appearing in our expansions are both homogeneous 
of degree —3 and smooth outside the origin. Therefore, our asymptotic profiles imply that 
conclusion (5) of Theorem 1 can be improved into 

U(x) = a(x) + 0(|x| -3 log \x\) , as|x|^oo, ifd = 2, 

and 

U(x) = a(x) + 0(\x\- 3 ), as|x|^oo, if d > 3 . 

The datum a can be replaced here by its filtered version e A a. 

It turns out that such improved estimates are optimal for generic self-similar solutions. 
For example, in the two dimensional case, the logarithmic factor cannot be removed, since the 
improved bound U(x) = a(x) + C(|x| -3 ) would require Q: A = 0: such stringent condition 
can be proved to be equivalent to the orthogonality relations J gl a\ = J gl a\ and J §i a\a2 = 0. 
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2.3 Main methods 



We will use the semigroup method and the theory of mild solutions of the Navier-Stokes 
equations as explained in detail in the books [4] and [13]. The main novelty of our approach 
relies on the use of the following ingredients: 

1. The first one is the use of remarkable, but not so much known, cancellation properties 
hidden inside the kernel K(x,t) of the Oseen operator e* A P, and inside other related 
operators, appearing in the integral formulation of Navier-Stokes. 

To be more precise, we can write K(x,t) = K(x) + t~ d / 2 K.2(x/\/t) , where &(x) is a 
tensor whose components are homogeneous function of degree —d (namely, second order 
derivatives of the fundamental solution of the Laplacian in M. d ), and K2 is exponentially 
decaying as \x\ — > 00. Such decomposition already played an important role in our 
previous work [3], where we showed that solutions u(x,t) arising from well-localized 
data behave like 



A crucial fact in the proof of the results of the present paper will be the use of the 
identities, for j = 1, . . . , d, 



Such cancellations are somehow hidden in K, because the non-homogeneous part K2 
(and, a fortiori, the kernel K) does not have a vanishing integral on the sphere. 

2. Our second ingredient are asymptotic formulae for convolution integrals: roughly speak- 
ing, these formulae consist in deducing the exact profile as \x\ -^00 of a convolution 
product f*g(x), from information on the regularity, the cancellations, and the behavior 
at infinity of the two factors / and g. In their simplest form, and for / and g "well 
behaved" at infinity, those formulae read 



We will apply several generalizations and variants of (7) in different situations (including 
the case of non-integrable functions) the factors / and g being either the Oseen kernel, 
the heat kernel, or a function related to the non-linearity. The assumptions for the 
validity of (7) are quite stringent (notice that (7) is obviously wrong if, e.g., f and g 
are both a gaussian function). Neverthless, the method that we use here has a wide 
applicability and can be used for constructing the far-field asymptotics for equations of 
other equations. See, e.g., [2] for an application to a class of convection equations with 
anomalous diffusion. 

We will also make use of the so called bi-integral formula. Such formula is obtained by 
simply iterating the usual integral formulation of the Navier-Stokes equations, which we now 
recall: 



Using the Oseen kernel K(x, t) , we can define the Navier-Stokes bilinear operator as 



u(x,t) ~ V x &(x) : E(t), as \x\ — > 00, 
where E{t) is the energy matrix of the flow: E{t) = (/ */ u^u^y , s) dy ds} . 





as \x\ — > 00 . 



(7) 




(NS) 



B{u,v){t)= / K(t- s) * V ■ (it® v){s)ds. 
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Then (NS) can be written simply as u = e a — B(u, u) . The bi-integral formula is obtained 
by a straightforward iteration: 

u (t) = e tA a - B{e tA a, e tA a) + 2B(e tA a, B(u, u)) - B(B(u, u),B(u, u)). (8) 

Roughly speaking, combining equation (8) with some nice properties of the heat kernel 
and fine decay estimates of the bilinear operator we can prove (e.g. when d = 3) that 
e tA a ~ a + tAa and B(e tA a,e tA a)(x) ~ tPV • (a <8> a)(x) as \x\ — > oo. After obtaining 
an explicit far-field asymptotics for u(x, i) , it is easy to deduce, in the self-similar case, the 
behavior at infinity of the profile U (x) , by passing to self-similar variables and eliminating t . 
The two last terms in the above bi-integral formula will contribute to the remaining terms in 
the right-hand side of expansion (6b). 

Notice that, in the two-dimensional case, the term PV-(a£3a) is not well-defined when a is 
homogeneous of degree — 1 . This explains the different structure of our asymptotic expansions 
in this case. The special structure of the asymptotic profiles in the two-dimensional case can 
be observed also if, instead of considering the behavior for \x\ — > oo as we do in this paper, 
one focuses on the the behavior of solutions for large time. (See, e.g., [8]). 

For sake of simplicity, in this paper we consider only data such that a(x) ~ M" 1 as 
\x\ — > oo, which is the natural assumption for the study of global strong solutions and related 
self-similarity phenomena. 

However, the study of the asymptotic behavior for large \x\ of solutions u (possibly defined 
only locally in time) is also of interest in more general situations, such as a(x) ~ |x| _1? - The 
far-field behavior of the solution u(x,t) of (NS), then mainly depends on the competition 
between three factors. The first one is the spatial localization of the datum (say, the value of 
the exponent $) and the consequent space-time decay of the linear evolution e tA a. The other 
two factors are the action of the quadratic non-linearity u(&u and of the non-local operator 
Pdiv(-) . 

When -d > d+1, the action of this nonlocal operator (whose kernel behaves at infinity like 
| x |-<i-i -j - g predominant, and is responsible of spatial spreading effects. When (d+ 1)/2 < $ < 
d+1 (the limit case •& = (d + l)/2 corresponding to the situations in which u®u decays like 
the kernel), the linear evolution becomes predominant and the spatial spreading phenomenon 
is not directly observed on the solution, but rather on its fluctuation u — e tA a. We refer to 
our previous paper [3] for a sharp description of these issues. 

The asymptotic profiles of u as \x\ — > oo in the cases < i? < 1 and 1 < i? < (d + l)/2 
should have a slightly different structure, but they are not known with precision yet. The 
method that we use in this paper for •& = 1 , and in particular the idea of iterating the 
Duhamel formula making use of the cancellations of the kernels, might be used to compute 
them. More and more iterations would be needed to deal with data decaying slower than 
| a? | 1 , or to determine the asymptotics to a higher order. On the other hand, no iteration or 
cancellation property was needed for the faster decaying data studied in [3]. 

The plan of the paper is the following: we begin with the study of the Oseen kernel. 
In Section 4, after some generalities about the asymptotic of convolutions, we describe the 
behavior at large distances of solutions to the heat equation. Section 5 is devoted to the (more 
or less standard) construction of solutions with a prescribed space-time decay. In Section 6 
we show how to use the cancellations of the Oseen kernel to get some new fine estimates. In 
the remaining part of the paper we will state and prove a more general form of Theorem 2. 



6 



3 Asymptotics and cancellations of the Oseen kernel DC and of 
the kernel F 

Let K(x,t) be the kernel of e tA P, let F{x,t) be the kernel of e* A Pdiv(-). Both K(-,t) and 
F(-,t) belong to C°°(R d ) and they satisfy the scaling properties K(x,t) = t~ d / 2 K(x/Vi,l) 
and F(x,t) = t-V+WFfr/y/i, 1) . 

Denote by T the Euler Gamma function and by 5j jk the Kronecker symbol. The following 
Proposition extends and completes a Lemma contained in [3]. 

Proposition 1 Let &= where &j tk (x) is the homogeneous function of degree —d 

6 M-M (S jtk \x\ 2 + dxjXk) 

KjM- 1 ) - 27r d/ 2 " \ x \ d + 2 ' y yd ) 

and J = (-Sj-h,k) j where dj-,h,k = dh&j,k, which we can write also as 

v M - r (^) a hh,k{ x )\ x \ 2 ~ {d + 2)x jXh x k 

#j;h,k(. X ) - " \x\ d+4 ' { ' 

where (Jj t h,k{x) = 5j,h,Xk + ^h,kXj + S k jXh , for j,h,k = 1, . . . , d. Then the following decom- 
positions hold: 

K(x,t) = St(x) + \x\~ d ^(x/y/tj, (10a) 

and 

F(x,t) =$(x) + \x\- d - 1 y(x/VtJ, (10b) 

where ^ and VP are smooth outside the origin and such that, for all a 6 N d , and x / 0, 
|3 a ^(x)| + |<9 Q, f(x)| < Ce -0 ^ . Here C and c are positive constant, depending on \ot\ but 
not on x . 

Moreover, the following cancellations hold: 

/ &(uj) da(u) = / uj£&(lo) da(uj) = 

J§<1-1 JSd-1 

[ ${u)da{uj)= I uj^{u)da{u) = (11) 
/ uj(U! m $(ui)d(j(ui) = 0, £,m=l,...,d. 

Remark 1 The homogeneous polynomials Q(x) = (Qj-h,k{x)) appearing in the statement 
of Theorem 2 is defined by the relation $(x) = |x|"^ 4 Q(x) , that is, with ^ d = T(^)/ir d/2 , 

Qj;h,k(x) = 7d(($j,hXk + S h;k Xj +S kJ x h )\x\ 2 - (d + 2)x j x h x k j. (12) 
Proof. The symbol of K is 

Taking the inverse Fourier transform we get 



djd k g s (x) ds = K^l(x, t) + K^(x, t) 



Computing the derivatives djd k g s (x) and changing the variable A = ^= in the integral we 
get 
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Kfl=7T- d / 2 \ X \- d 
j,k l l 



\x\/VU 



J ' \x\ 2 



dX. 



But, for all r > and a > — 1, 
Choosing first a = d — 1 , then a = d + 1 and using T((d + 2)/2) = (d/2)r(<i/2) , we get 



X 

2^ a 



Oj,k + » , lo 



+ \x\- d V j:k (x/V4t). 



Here, ^ = (^j.fc) is a family of functions such that, 

VaGN d , \d a ^{y)\ <C a e~ c ^ 2 , y G R 3 . 



(13) 



Observing that IK^ can be bounded by the second term on the right hand side and modifying, 
if necessary, the functions ^j,k (which can be done without affecting estimate (13)) we see 
that decomposition (10a) holds. The decomposition (10b) is now an immediate consequence 
of the definition of $(x) . 

Observe that = djO^Ed, where E^ is the fundamental solution of —A in W 1 . From 
the radial symmetry of E^ , we immediately get 

/ &j k(uj) da(u>) = / u>j&(uj) da(uj) = 0, j / k 

J§d-1 ' J§ d -! 



and 



/ $(w)da(u) = / weuj m $(uj)da(uj) = 0, 

Jgd-l Jgd-1 



£, m = 1, . . . , d. 



Using again the radiality of E^ and AE^ = on S d 1 , yields f S d-i da{io) = 0. This 

argument also shows that the identities 



/ 

J§d-1 



ui $j-h,k{u) da(u) =0, j,h,k,£ = l,...,d 



can be reduced to the proof of the equality 

/ Lu i d e d]E d (u)da(Lu)= [ u £ d$E d (cj)da(u;), j^L (14) 

The fact that both terms in (14) are zero follows from djdhdkE^uj) = Qj i / l ,fc(<^) , for uj 6 
and formula (12). In the computation, one needs to use the moment relation 



/ uPl d(j{uS) = — f da{uj) 

JSd- 1 d J§d-l 



and the well known identities (easily obtained via the Stokes formula) 

3 



cjj da {to ) 



dcr(uj) 



d(d + 2) J S d-i 

f uj]ulda(uj) = ,,,\ , f da(uj), j / k. 
I J§d-i d{d + 2) J S d-i 



□ 
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4 Far-field asymptotics of convolutions and application to the 
heat equation 

The purpose of our next result is to describe the exact behavior as \x\ — > oo of the convolution 
product of two functions / and g from the asymptotic properties of each factor. We will 
consider only a simple particular situation that will be sufficient for our purposes. 

Proposition 2 Let d > 1 and rn > 1 two integers. Let f G E™ for some < $ < d and 
g G L 1 (M d , (1 + \x\) m dx) n -E^+m • Then the convolution product f * g satisfies 



f*g(x) = ^f~(/ y 7 5(y)rfy)^/(x) + n{x), 



(15) 



7£ 

0<M<m-l 



where 1Z{x) is a remainder term satisfying, for some constant C > independent of f and 
g and all x / : 

\n(x)\ < ckl" m " ,? ll/ll^(lbbo +m + \\g\\ L i(R*,\ x \™*c))- (is) 

Remark 2 The identity (15) is useful, for large \x\ , when at least one derivative |3 7 /| decays 
at infinity exactly | 7 (at least in some directions). In this case, lZ(x) is indeed a 

lower order term as |x| — > oo. 

Proof. We can assume, without restriction, that ||/||pm = \\g\\po = 1- We have to estimate 

the difference between J f(x — y)g(y)dy and the first term on the right-hand side of (15). 
Such difference can be written as the sum of four terms D\ + ■ ■ ■ + D4 , where 

Di= [ \f(x-y)- Yl { - Z ^dy(x)y^}g(y)dy, 

D 2 = g(x -y)f(y)dy, 

J\v\<\*\/2 

D 3 = I f(x- y)g{y)dydy 

J\y\>\x\/2, \x-y\>\x\/2 

and 

Di = _ ^j^dy(x) [ y^g(y)dy. 

i 7 i<m-i 71 J\v\>\*\/2 

Using the Taylor formula, we see that 

|A|<C|x|-*- TO / \y\ m \g(y)\dy, 

which is bounded by the right-hand side of (16). Direct estimates show that I-D2I, \D$\ and 
I-D4I are bounded by C\x\~®~ m as well. 

□ 



Remark 3 We can give now a more precise statement about the asymptotics claimed in (7). 
The simplest result reads as follow: if f,g £ E^ +d (the non-homogeneous space) for some 
a > 0, then 

/ * g(x) = ( J f)g(x) + ( J g)f(x) + 0(\x\- d - a *), 
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as \x\ — > oo, where a* = min{2a, a + 1}. When a = 1 the remainder must be replaced 
by (D(\x\~ d ~ 2 log(|x|)) . The proof relies on the same argument as that used in the proof 
of Proposition 2: the only difference is that the Taylor formula is applied to both / and 
g , so that one has to introduce an additional term D§ in the decomposition of / * g . Of 
course, one could state many variants of this result: here the most important condition was 
that the decay of the two factors / and g (or at least the decay of the factor with the least 
spatial localization) must increase after derivation; but one could put, instead, a more general 
condition in terms of moduli of continuity. 

Other useful functional spaces are, for m £ N, i? > 0, 
Xfi 1 = (tie Lj oc ((0,oo),C m (R d )): \\u\\ x ? = max esssup (V~t + |x|f +H \d%u{x, t)\ < oo). 

I ^ \a\<m x ,t J 

The use of such spaces for the Navier-Stokes equations is more or less classical (see, e.g., 
[6], [7]) but, unfortunately, there is no agreement on the notations. 

The following lemma is elementary: 

Lemma 1 Let m £ N, a £ E™ , with < d < d. Then there is a constant C > 0, 
independent on a, such that 

ll e * Aa llxT < C||o||pm. 

Proof. For a G N d , \a\ < m, one writes d^e tA a(x) = J d%gt(x — y)a(y)dy and splits 
the integral in M. d into the three new integrals, corresponding to the three disjoint regions 
\y\ < M/2, [a? — 2/| < \x\/2 and the complementary region in M. d . For the second integral one 
first applies \a\ -times integration by parts. Then the direct estimate \d%gt(x)\ < C\x\~ d ~\ a \ 
gives the spatial decay |9"e iA a(x)| < C | a; | ^ I "I . On the other hand, a belongs to the Lorentz 
space L d /*>°°(R d ), and d*g t G L d /( d -*)' 1 (R d ) . Then the time decay estimate ||^e tA a||oo < 
Qj.-(ti+\a\)/2 follows from the generalized Young inequality (see, e.g. [13]). 

□ 

As an application, we get the exact asymptotic profile as \x\ — > oo for the solution of the 
Cauchy problem associated with the heat equation for slowly oscillating data. We first recall 
two standard notations: if (3 £ N d we set: (2(3 - 1)!! = Y[j=i,...,d 1 • 3 • . . . • (2(3j - 1) and 

/3,->l 

(2/3)!! = Y\j=i,...,d 2 ■ 4 • . . . • 2(3. Now we can state the following: 

0j >i 

Lemma 2 (i) Let m > 1 he an integer, < i? < d and a £ E™ . Then, 

e tA a(x) = Yl ^m\ dWa ( x) + (* m/2 M _,? ~ m )> as N oo, 

|2/3|<m-l ^ ' '" 

uniformly for t > (i.e. the remainder term is bounded by Ct m l 2 \x\~' d ~ m ). 
(ii) In particular, if m > 4 and a £ E™ : 

e tA a(x) =a(x) +tAa(x) + C(t 2 |x|" 5 ), as \x\ -> oo, 

uniformly for t > . 

Proof. Indeed, writing e tA a(x) = g t * a(x) , we can apply Proposition 2 with gt(x) = 
(Ant)- d / 2 e-\ x \ 2 /^ instead of g. Observing that, for all (3 £ N d , 

J y^g t (y)dy = 2^(2P-l)\\t^ 

and that f y 1 gt(y) dy = if 7 £ N d is not of the form 7 = 2(3 , we obtain the result. 

□ 
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5 Global existence of decaying solutions 

We already recalled that F denotes the kernel of e* A Pdiv(-) and, for t > 0, F(x,t) = 
t-( d+1 V 2 F(x/Vi,l). It is also well known that \d£F(x,l)\ < C a {l + | £c | ) _cf_ 1_ for all 
a G N d . A quick way to prove this decay at infinity is to observe that such estimate is 
immediate for \x\ > 1 for both terms in the right-hand side of equation (10b). Moreover, it 
is clear from its definition that F(-,t) £ C°°(R d ) for t > 0. 

Let us introduce the linear operator L, defined on d x d matrices w = (wh t k) by the 
relation ^ 

L,(w)(x,t) = J J F(x - y,t - s)w(y,s)dyds. (17) 

More explicitly (and accordingly with the notation introduced in Section 2), the j -component 
is given by 

h(w)j(x,t) = / y^Fj-h^x - y,t - s)w h)k {y,s)dyds. 
JoJ h.k 

The interest of considering such operator is that the Navier-Stokes bilinear operator can 
be expressed as 

B(u,v) = h(u®v). 

We start with a simple lemma (already known in a slightly less general form, see [15], [7]). 

Lemma 3 Let m £ N and w = (wh : k) £ • Then h(w) £ X™ and, for some constant 
C > independent of w, 

\\Uw)\\x r < C\\w\\ XT . (18) 
Proof. We can assume, with no loss of generality, Hu'llxj 1 = 1- We start writing 

dZL(w)(t)=J*J b^F(x-y,t-s)w(y, S )ds. (19) 

Let a G N d , such that |a| < m and x / 0. We split the spatial integral in equation (19) into 
the three regions |y| < \x\/2 , next \x — y\ < \x\/2, then ( \y\ > \x\/2 and \x — y\ > \x\/2) and 
we denote with I\ , I2 and ^3 , the three corresponding integrals. Prom the estimate (deduced 
from (10b)) \d%F{x,t)\ < C\x\-( d+ \ a \h~ 1 / 2 , and the estimate \w(y,s)\ < ^^s- 1 / 2 , we 
obtain immediately 

\h(x,t)\ + \h(x,t)\<C\x\~ l -\ a \. (20) 

We now treat I2 ■ When \a\ = we can simply use well known fact that ||F(-,t)||i < Ct~ x l 2 
to obtain \l2(x,t)\ < Cjxl™ 1 . When 1 < |a| < m we make as many integration by parts as 
needed, and use estimates of the form (deduced from the rescaling properties of F recalled 
at the beginning of this section and the fact that d%F(-, 1) G L 1 (M. d , (1 + |x|)l a l dx) ) 

11 1 - r^^c-, * - S ) iii < c(* - S )-V2. 

Then observing that \dyiu(y, s)\ < \y\~ 1 ~^s~ 1 ^ 2 for \a\ < m, we conclude that I2 can be 
estimated like I\ and ^3 in (20). Summarizing, we showed that 

|S£L(u;)(M)| < Clxr 1 " 1 " 1 . (21) 

There is a now well known strategy (see [15]) to deduce time decay estimates from the 
corresponding space decay estimates. Namely, using the semi-group property of the Oseen 
kernel, 

h(w)(t) = e tA/2 L(w)(t/2) + [ F(t-s)* w(s) ds = K±{t) + K 2 {t). 

Jt/2 
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From the Young inequality in Lorentz spaces, and observing that ||L(w;)(£)|| id:00 is uniformly 
bounded, because of inequality (21), we get 

||g?*i(t)L <||^/2|| J/V(i -i,,i||L(u;)(t/2)|| LAoo < ct-^\y\ 

Moreover, 

ll^^wiu < f \m-s)u\^w{8)Uds<ct-^\y 2 . 

Jt/2 

Concluding, we showed that 

\d%L(w)\(x,t) < C(\x\-( 1+ \ a V aH 1+ H)/ 2 ) < C'{Vt + \x\)- l -\ a \. 
This proves the natural estimate (18). 

□ 

We now follow the standard procedure for constructing global solutions to (NS) in the 
space X™ . Our starting point will be the following basic existence result, which is nothing 
but a reformulation of well-know results in the literature (see [4], [5], [7], [15]) in a slightly 
more general form. 

Proposition 3 Let d > 2 and m > he two integers. There exist two constants e > and 
M > such that for all divergence-free vector field a 6 E™ , satisfying 

&|| g™ ^ e > 

there exists a unique solution u G X™ of (NS) starting from a (in the sense that u(t) — > a 
in S'(M. d ), as t -> OJ, such that < eM. 

Proof. We only have to apply the size estimate for the linear evolution 

||e tA a||xf < C||a||^m 

(this is a particular case of Lemma 1) and the corresponding estimate for the bilinear operator: 

\\B(u,v)\\ x ™ < C||u||xj"||v||xf- 

This last inequality is obtained applying Lemma 3 with w = u (g> v . The existence a solu- 
tion u £ X[ n (and its unicity in a ball of such space) now follows from the application of 
the contraction mapping theorem, as explained e.g. in Cannone's book [4]. Slightly chang- 
ing the estimates of the previous Lemma we easily obtain, e.g., the bound \B(u,u)(x,t)\ < 
C\x\^l 2 t 1 ^ , implying B(u,u)(t) -> in S'(R d ) as t -> 0. Thus, from (NS), u(t) -> a as 
t — > in the distributional sense. 

□ 

Remark 4 In the particular case in which a is a homogeneous vector field of degree — 1 
in M. d , the solution u constructed in Proposition 3 is self-similar: 

U(M) = ^(^)' 

for some with U G E 1 ™ (the non-homogeneous space). This easily follows from the scaling 
invariance of (NS) (see e.g. [4, Ch. 3]). 
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6 Fine estimates of the bilinear term 

It follows from Lemma 3 that, for w G X®, we have 

|L(w)(x, t)\ < C{r 1 / 2 A ^p 1 ). (22) 

This was enough for constructing a decaying solution of (NS). 

However, to obtain such decay estimate we used only few properties of the kernel F(x, t) , 
namely, its pointwise decay and its rescaling properties. Next Lemma will allow us to consid- 
erably improve estimate (22), at least in the parabolic region \x\ > y/i. Its proof will make an 
essential use of the cancellations properties of the kernel F(x, t) and requires some regularity 
for w . 

Lemma 4 Let w = (wh : k), with w G X 2 ■ Let h(w) be defined by equality (17). Then we 
have, for d > 3 , 

\h(w)(x,t)\ <c(t- l/2 A t\x\~^. (23a) 
When d = 2, we have the weaker estimate 

\h(w)(x,t)\ < Ct|x|~ 3 log(^i), \x\>eVi. (23b) 

Under the more stringent assumption w 6 Xf , we have the following estimates for VL(w) : 

'C{t- 1 At|x|~ 4 ), ifd>3 
|VL(i«)(x,t)| < < Ct- 1 ifd = 2 and \x\ < e^ft 

Ct\x\- A \og(\x\/^ft) ifd = 2 and \x\ > e^ft. 

In all these inequalities C > is a constant dependent on w only through its || ■ \\ X 2 or its 
|| • || x z -norm, and independent on x and t. 

Proof. We can limit ourselves to the region \x\ > e\ft. Indeed, when \x\ < ey/i the result holds 
because of inequality (18), which, in the special case m = 0, 1, implies |L(u;)(x, t)| < Ct~ x l 2 
and |VL(x)(x,t)| < Ct' 1 . 
Let us decompose 

h(w)(x,t) = / / F(x — y,t — s)w(y, s) dy ds 

JoJ\y\<\x\/2 

+ ii F(y,t- s)w(x-y,s)dyds 

JoJ\y\<\x\/2 (24) 

+ // F(x-y,t- s)w(y,s)dyds 

JoJ\y\>\x\/2, \x-y\>\x\/2 

= Li + L 2 + L 3 

We start with estimating L3. Using \F(x — y,t — s)\ < C\x — y^^ 1 < C'\y\~ d ~ l (the two 
inequalities being valid in the region of M d where we perform the integration) and \w(y, s)\ < 
\y\~ 2 , we get |L 3 (x,t)| < Ct|x|" 3 . 

In view of the use of the Taylor formula, we further decompose Li (recalling also (10b)) 

as 

Li = / / [F(x — y, t — s) — F(x, t — s)] w(y, s) dy ds 

JoJ\y\<\x\/2 

+ $(x): [ [ w(y,s)dyds (25) 

Jo J\y\<\x\/2 
ct 



+ \x\~ d 1 [ ^(x/Vt - s) [ w(y,s)dyds. 

Jo J\y\<\x\/2 
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Using \VF(x,t)\ < C\x\~ d ~ 2 , next \y\ \w(y,s)\ < C|y| _1 shows that the first term in (25) is 
bounded by Ct|x|~ 3 . 

When d > 3, since \w(y, s)\ < C\y\~ 2 , the second term in the right-hand side of (25) is 
also bounded by Ct\x\~ 3 . When ci = 2,make use of the inequality \w(y,s)\ < C + \y\)~ 2 
and of the chage of variables y = \fsz . This leads to the weaker upper bound estimate of the 
form Ct\x\~ 3 log(\x\/y/i) , valid for \x\ > ey/i. 

The simplest way to treat the third term on the right-hand side of (25) is to recall that 
|^(x)| < C . In this way, one can proceed exactly as for the previous term and obtain the 
same bounds. This would be enough for the proof of this Lemma. However, for later use 
(namely, to shorten the proof of Lemma 6 below), we want to prove that this last term in (25) 
is bounded, in the region \x\ > ey/i, by Ct|x| -3 also when d = 2. This is easy: indeed ^ has 
a fast decay at infinity; here, the use of the inequality |*(x)| < Cjxj" 1 is enough to conclude. 

We now consider L2 . We decompose it as 

L 2 = / / F(y,t — s)[w(x — y, s) — w(x, s) + y ■ Vw(x, s)] dy ds 

JoJ\y\<\x\/2 

ft 



[ w(x,s) [ F(y,t-s)dyds (26) 

Jo ^|y|<|x|/2 

-/ Vw(x,s)- yF(y,t- s)dyds. 

Jo J\y\<\x\/2 

Now we use the inequalities \V 2 w(x,t)\ < C\x\~ 4 and \y\ 2 \F(y,t — s)\ < C\y\~ d+1 , and 
obtain that the first term on the right-hand side in (26) is bounded by C£|x| -3 . We now 
conclude using the cancellations of the kernel F: more precisely, since J F(-,t — s) dy = 
and \F(y, t - s)\ < C^-^ 1 the second term is also bounded by Ct\x\ . 

A brutal estimate of the third term in (26) would give a non-optimal bound of the form 
C|x|~ 3 log(|x|\/i) for large \x\, which is not enough. But, for \x\ > 2y/i, the third term 
in (26) can be further decomposed as 



/ Vw(x,s)- _yF(y,t- s) dy ds 

Jo Jlvl^V^s 

+ f Vw{x,s)- f y$(y)dyds (27) 

Jo JVt :r s<\y\<\x\/2 

+ / Vw(x,s)- [ y\y\~ d ~ 1 ^(y/y/t zr s)dyds. 

Jo Jx/t^<\v\<\x\/2 



Now it is easy to see that the first and the third term in (27) are 0(t|x|~ 3 ). But 5 has 
vanishing first order moments on the sphere (see Proposition 1) so that the second term 
in (27) is zero. 

Summarizing, we have established inequality (23b) in the two-dimensional case and in- 
equality (23a) when d > 3 . 

To prove the inequality for VL, we fix I G {1, . . . , d} and we write 

diL(x, t) = / F(x — y,t — s)d(w{y, s) dy ds = Li + L 2 + L 3 , 
Jo 

where the decoposition is obtained as before (see (24)). The two terms L 2 and L3 are treated 
exactly as before, but we get now upper bound of the form Ct|x|~ 4 since dew (and its 
derivatives up to the second order) decays faster than w (and its corresponding derivatives). 
Notice that we need use here the assumption w 6 Xf which ensures a decay for the derivatives 
up to the order three. 
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For treating Li we integrate by parts. It is easy to see that the boundary term is bounded 
by Ct\x\~ 4 . The other term is J* * f\ y \<\ x \/2 deF(x — y, t — s)w(y, s) dy ds , for which we obtain 
the usual bound Ct|x| -4 when d > 3 and Ci|x|~ 4 \og(\x\/y/i) for d = 2 and \x\ > e\ft. 

□ 

Remark 5 For later use, let us observe that if u € X 2 is the solution constructed in Propo- 
sition 3, in the case m > 2, then, applying Lemma 4 to w = u (g> u, so that L(to) = B(u, u) , 
we get 

'C(t-V2 Ai|x|~ 3 ) if > 3 

Ct -1 / 2 if = 2 and |x| < (28a) 

Ct\x\~ 3 \og{\x\/^) if d = 2 and |x| > e>/t. 



In the case it G Xf (this requires the more stringent assumption a £ Ef in Proposition 3), 
in addition to the above estimates, the bilinear term satisfies 



\VB(u,u)\(x,t) < < 



'CQ- 1 At\x\~ 4 ), if d > 3 

Ct~ l lid = 2 and \x\ < e^t (28b) 

i Ct|a;|- 4 log(|x|/Vt) if d = 2 and |x| > eVt. 



These estimates will play an essential role in the study of the bi-integral formula 

u (t) = e tA a - B(e tA a, e tA a) + 2B{e tA a, B(u, u)) - B(B(u, u),B(u, u)). (29) 

7 Asymptotic profiles of the velocity field in the 2D case 

In the two-dimensional case, from Lemma 1 and Remark 5 we get, for (x,t) £l 2 x (0,oo), 
, fA x / m I Ct^ 1 if Ixl < e\/t . , 

i e - o0B(u , .)(,, t) i s { ctW - llog(|x|M) if ;j; eVi . (so) 

The last term in (29) satisfies, always for (x, i) G M 2 x (0, oo) , an even stronger estimate, 
namely 

it,/ m fc* _1 if M < e\/t 

i«(».»)«^.»)(->*)i<{ ctV| - 6log2(k|M) ,^1;^, mi 

Next Lemma allows us to show that the two last terms in the right-hand side of (29) can 
be considered as remainders, i.e., they can be included in the C(t|x|~ 3 ) term. 

Lemma 5 Let w = (wh^k) defined on M 2 x (0, oo) with Wh,k(x, t) bounded by the right hand 
side of (30), or by the right-hand side of (31). Then, if h(w) is given by (17), we have for 
some C > independent on x or t, 



\Ux,t)\ < eft' 1 / 2 A i|x|- 3 ). 



Proof. Our assumptions imply w £ X®. Then we deduce from Lemma 3 that |L(x,t)| < 
Cr 1 / 2 , therefore we can assume that |x| > ey/t. Then we split the spatial integral defining 
L (see (17)) into the three regions \y\ < y/s, \fs < \y\ < |x|/2 and |y| > |x|/2. The first 
term that we obtain is bounded using \F(x — y,t — s)\ < C|x| -3 (this is true only in 2D) 
and \w(y, s)\ < Cs -1 . For the second term we use the same bound for F and \w(y, s)\ < 
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Cs|y| 4 log(|y|/y / s) . The last term is treated using the bound \w(y,s)\ < C\fs\y\ 3 and that 
mt - s ) h <C{t-s)-V\ 

□ 

Next Lemma will be useful for treating the term B(e tA a, e tA a) arising in (29). Note that 
for a G E\ we have, from Lemma 1, e tA a ® e* A a G X| . 

Lemma 6 Let w = (wh : k) , with Wh : k £ X% for all h,k = 1,2. Then we have 

L(w)(x,t) =$(x): [ [ w(y,s)dyds + 0(t\x\~ 3 ), as \x\ -> oo, (32) 

uniformly with respect to t in the region \x\ > ey/t . Here h(w) is given by (17) and $(x) is 
the homomeneous tensor of order three defined by equation (9b) . 

Proof. This follows from the proof of Lemma 4. Therein, we decomposed h(w) as the sum of 
several terms, all of which, excepted one, could be bounded by Ct|x| -3 . The only term for 
which such upper bound could brake down was 



$(x): [ [ w{y,s)dyds 

JO J\y\<\x\/2 



10 J\y\<\x\/2 

(see the second term in the right-hand side of (25)). A simple modification of the error term 
now shows that we can change the above domain of the spatial integral into {\y\ < \x\} . 

□ 

Lemma 7 Let a(x) be a vector field defined on M? , such that a G E\ . Then, for \x\ — > oo 
and uniformly in time, in the region \x\ > e\ft, we have: 

rt 

sA„ ^ „sA 



{e s£X a®e s£X o)(y)dyds= / (a (g> a){y) dy ds + 0(t 1) 

0^|»|<|x| Jo J^<\y\<\x\ 

(here and below 0(tl) denotes a remainder function bounded by Ct for \x\ > ey/i). 
In particular, if a is homogeneous in M 2 of degree — 1 



f f {e sA a®e sA a)(y)dyds = tlog(^)( f a®a)+0(tl), 
JoJ\y\<\x\ V V* y KJS 1 J 



as \x\ — > oo . 



Proof. Indeed, we can assume \x\ > Then / / (e s a®e s a)(y)dyds is bounded 

JoJ\y\<^ 

by Ct. It remains to treat 

/ / (e sA a®e sA a)(y)dyds, 
Jo J^<\y\<\x\ 

which we can rewrite as the sum of four new integrals, if we use the decomposition e tA a(x) = 
a(x) + 1Z(x,t) obtained in Lemma 2 (in the case $ = 1, m = 2) and a similar decomposition 
for e tA b. Here, K satisfies \K(x,t)\ < Ct\x\~ 3 . An easy calculation shows that the three 
integrals containing at least one factor 1Z are bounded by Ct. 

□ 
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Theorem 3 Let u(x,t) £ X\ he the global solution of the Navier-Stokes equations in M 2 , 
with datum a € E\ (as constructed in Proposition 3). Then u has the following profile for 
\x\ — ► oo, uniformly with respect to t in the region \x\ 

u(x, t) = a(x) - $(x) : [ [ (a ® a)(y) dy ds + 0(t\x\~ 3 ). (33) 

Moreover, if a is homogeneous of degree — 1, then u(x,t) = -^U^-^j is self-similar and the 
profile U (x) is such that 



U(x) = a(x) — log(\x\)$(x) : 



(^J a (g> aj + 0(\x\ 3 ), as|x|^oo. (34) 



Proof. The first statement follows from the bi-integral formula (29) and our previous Lem- 
mata. Indeed, as we have already observed, by Lemma 2, we can write e tA a(x) = a(x) + 
C(t|x|" 3 ). Next, writing 

B(e tA a,e tA a) = Ue tA a®e tA a), 

we apply first Lemma 6 with w = e tA a ® e tA a, and then Lemma 7. This shows that 
—B(e tA a,e tA a) equals to the second term on the right hand side of (33), up to an error 
0(t|x| -3 ) for large \x\ . The last two terms in the bi-integral formula can also be included into 
the remainder term 0(i|x| -3 ), as shown by combining inequalities (30)-(31) with Lemma 5. 
In the case of homogeneous data, an elementary computations shows that 

/ / (a®a)(y)dyds=([ a <g> a) t logf ^) + t/2. 

JO J \/s<\y\<\x\ KJS 1 J V V* y 

Then profile (34) follows from profile (33) passing to self-similar variables and eliminating t . 

□ 

8 Asymptotics in the higher-dimensional case 

We now establish the analogue of Lemma 6 for the higher dimensional case. 

Lemma 8 Let w = (wh,k) with Wh,k £ X\. Then we have, as \x\ —> oo, and uniformly in 
time, for \x\ > e\ft , 

L(w)(x,t) =$(x):J^J w{y,s)dyds + 0(t\x\~ 5 log(\x\/Vt)) (35a) 

for d = 3, and 

L(w)(x,t) = 0(t\x\~ 5 log(\x\/Vi)) (35b) 

when d > 4 . 

Proof. We go back to the decomposition L = Li + L2 + L3 obtained in (24). Writing Li as 
in (25) and using the estimate \w(y,s)\ < C(y / s + |x|) -4 , the bound \VF(x,t)\ < C\x\~ d ~ 2 , 
and the fast decay of \I> shows that the first and the third term in (25) are bounded by Ct|rc|~ 5 
(with an additional logarithmic factor log(|x|/\/t) , for the first term in (25), when d = 3) for 
\x\ > e\/t. The second term in (25) has the form 



m-- ff 

JoJ\y 



w(y, s) dy ds. 

l<N/2 
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Using again that \w(y,s)\ < C{y/s + | a; | ) 4 and distinguishing the between the cases d = 3 
and d > 4 shows that such term can be written as the right-hand sides in (35a)-(35b). 
We now decompose L2 , as 



L 2 = / / F(y,t — s)[w(x — y,s) — w(x,s)] dyds 

JoJ\y\<\x\/2 

+ / w(x,s) / F(y,t- s)dyds. 

Jo J\y\<\x\/2 



'\y\<\x\/2 

Since \Vw(x,t)\ < C|x| -5 , the first term in (36) is bounded by C\x\~ 5 tlog(\x\/\/t) for \x\ > 
e\ft. Combining the estimate \F(y,t — s)\ < 1 2/ 1 rf 1 with the condition f F(-,t — s) ds = 0, 
shows that the second term in (36) is bounded by C£|x| — 5 . Such bound holds also for L3 as 
easily checked using the usual spatial decay estimates of F and w . 

□ 

Our next Lemma essentially states that if a and b are two functions defined on M. d and 
well behaved at infinity (for example, the derivatives of a and b decay faster than a and b 
as \x\ — > 00), then 

(e tA a)(e tA b) ~ e* A (a&), as \x\ -► 00. 

More precisely, we have: 

Lemma 9 Let d > 3 and a,b € E\ . Then 



{e tA a){e tA b) = e tA (ab) - 2 f e^ 3 ^ [Ve sA a • Ve sA b] ds. 

Jo 



(37) 



Proof. Let v = e tA a and w = e tA b. Then we have dtv = Av and dtw = Aw. Multiplying 
by w the first equation and by v the second one we get 

dt{vw) = wAv + vAw = A(vw) — 2Vv • Vw. 

Since d > 3, ab is locally integrable in M. d . But (vw)(t) — > a6 as t — > + weakly (because 
■u(i) — > a and ■«;(£) — > 6 in L^ oc (M^) , for example, as t > 0). Then the conclusion follows 
from Duhamel formula. 

□ 

In the above Lemma we only used, in fact, a,b G E® . The stronger assumption a,b G E\, 
however, ensures that the last term in (37), decays faster as \x\ —> 00 than e tA (ab). 

We now give the higher-dimensional counterpart of Theorem 3. 

Theorem 4 Let u(x, t) G Xf be the global solution of the Navier-Stokes equations starting 
from a G Ef (as constructed in Proposition 3). Then u has the following profile as \x\ —> 00, 
uniformly in time for \x\ > c\ft . For d = 3, 

u(x,t) = e tA a(x) -te tA PV- (a ® a) -ff(x):A(t) + O ^ 2 |x|~ 5 log^)^ , (38a) 

for some matrix-valued function A(t) = (A^fe(i)), satisfying \A(t)\ < Ct 3 ^ 2 . Moreover, when 
d> 4, 



u (x,i) =e* A o(x) -te* A PV- (a^a) + 0^ 2 |xr 5 log(^=)^). 



(38b) 
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Remark 6 The function A(t) is not know explicitly, but it depends on u and a in an explicit 
way: see formula (41) below. For more regular data, namely a 6 Ef , and recalling Lemma 2 
(applied with m = 4 and ■d = 1) one can replace in the above asymptotics the term e tA a(x) 
with a(x) + tAa(x) . 

Proof. As for the proof of our previous theorem, we write u by means of the bi-integral 
formula (29). As an application of Lemma 9 we can rewrite (for d > 3) the term B(e tA a, e tA a) 
appearing in the bi-integral formula (29) in a more convenient form (we denote here by T A 
the transposed of the matrix A): 



(e tA a,e tA a)= [ F(t - s) * (e sA a <g> e sA a) ds 
Jo 

= J e^^ A Fe sA V ■ (a ® a) - 2 J J e^^PV • T (v ® e rA aj (y ® e rA a] 



te tA F ■ V(a ® a) - 2 (t - r)e (( - T)A PV • e rA a) (^V e rA a) 



dr , 

(39) 

where we applied Fubini's theorem in the last equality. 

We set 

L(w)(t)= [ (t — t )F(t — t) * w(t ) dr and h(w)(t) = [ r F(t - r) * w(t) dr. 
Jo Jo 

Note that, excepted for the additional factors t — r or r, the operator L and L agree with 
the operator L introduced in (17) and studied before. If we introduce the matrix 

wi = T (V ® e rA aj (V tg) e rA aj , 

then we can rewrite (39) as 

B(e tA a, e tA a) = t e tA FV ■ {a <g> a) - 2L(u>i). 

The estimates of Lemma 1 (in the case m = i? = 1), imply w\ G X| . But the result of 
Lemma 8, established before for the operator L, can be easily adapted to the operators L 
and L; indeed the factors t — r and r are harmless in our estimates due to the obvious 
inequalities t — r <t and r <t. Thus, we get, for d = 3, 

L(w 1 )=!5(x):J(t-T)Jw 1 dyds + O^Vr 5 log(|x|M)) . 
When d > 4, we can simply write 

L(^) = 0(t 2 |a:r 5 log(|x|/V£)). 

It remains to write the asymptotics (or to estimate) the two last terms B(e tA a, B(u,u)) 
and B(B(u,u), B(u,u)) appearing in the bi-integral formula (29). Let 

W2 = j e tA a (8> S(u, n). 

We get from Lemma 1 (applied with m = 1 and $ = 1 ) and Remark 5 that W2 £ - In the 
same way, Remark 5 ensures that, if we set 

ti>3 = j B(u, u) <g> B(u, u), 
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then W3 £ X\ . Therefore, Lemma 8 (or more precisely, the adaptation of this Lemma to 
h(u)2) and L(u>3) ) implies, for d = 3, 



2B(e sA a,B(u, u)) - B(B(u, u), B(u, u)) 
= 2L(w 2 ) - h(w 3 ) 

= $(x)--J sJ{2w 2 -w 3 )dyds + o(t 2 \x\- 5 log(\x\/Vt)^ 



as \x\ -^00. 



When d > 4, the first term in the right-hand side of (40) can be dropped. Therefore, 
the proof of the expansion (38b) follows from the bi-integral formula, collecting the above 
estimates. 

In the case d = 3 , it is convenient to introduce the time-dependent matrix 



A(t) = J J -2(t - s)wi - 2s w 2 + SW3 



dyds. (41) 



The expansion (38a) now follows by collecting all the above expressions. The estimate 
|A(i)| < Ct 3//2 is immediate, because wi , w 2 and W3 belong to X\ . 

□ 

As an application of this theorem, we can complete the proof of Theorem 2 by giving the 
far-field asymptotics of self-similar solutions in the case d > 3 . 

End of the Proof of Theorem 2. We assumed that a G C OD (§> d ~ 1 ) and that a is is homogeneous 
of degree —1. From the second part of Lemma 2, 

e tA a{x) = a(x) + tAa(x) + 0(t 2 \x\- 5 ). 

But the solution u is of the self-similar form u(x,t) = -^U(xj\fi). Moreover, the linear part 

e tA a and the nonlinear part B(u,u) of u are also of self-similar form, so that, with the same 
notations of the previous proof, Wj(y,s) = jpWj(y / \/s) , where 

W j (y)=w j (y,l), j = 1,2,3. 

If follows from (41) that, in the case d = 3, A(i) is of the form A(t) = t 3 / 2 B, for some 
constant matrix B = (Bh,k) ■ As for A(t), such matrix B is not known explicitly, however, it 
is possible to obtain an explicit integral formula relating B to the datum a and the profile 
U, performing a self-similar change of variables in the integral (41). An easy computation 
yields 

B = \f (-swi - + 2W z) (y) d v- C 42 ) 

Now we can pass to self-similar variables in expansion (38a) and, after eliminating t, we 
get, for d = 3, 

U(x) = a(x) + Aa(x) - e A P • V(a a)- +0(\x\~ 5 log(|x|)), (43a) 

as \x\ ^00. 

As before, for d > 4, the far-field asymptotics has a simpler structure, namely, 

U{x) = a(x) + Aa(x) - e A P • V(o ® a) + 0(\x\- 5 log(|x|)) , (43b) 
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as \x\ — > oo . 

To finish the proof, it remains to show that we can drop the filtering operator e appearing 
in the right-hand side of equations (43a) and (43b). Recall that a is smooth on the sphere. In 
fact, the condition a G C°°(S a!_1 ) will allow us to carry the proof using only "soft arguments". 
The datum a being homogeneous of degree — 1 , V ■ (a <8> a) is a homogeneous distribution of 
degree —3 (here we need d > 3), which agree with a C°° function outside the origin. But 
the matrix Fourier multiplier of the operator P (given by 5j : k — £j'6fc|£l _2 ) 1S homogeneous 
of degree zero and smooth outside the origin). Then it follows (see, e.g., [17, p. 262]) that 
PV ■ (a <8) a) is a homogeneous distribution of degree —3 that agrees with a C°° function 
outside the origin. 

Now let x ^ C^°(IR d ) be a cut-off function equal to 1 in a neighborhood of the origin and 
write 

e A PV ■ (a <g> a) = e A X PV • (a <g> a) + e A (l - x)A(x), 

where A(x) a smooth function on M. d , agreeing with PV • (a ® a) outside a neighborhood of 
the origin. In particular, (1 — %).4 G E™ , for all m G N. 

Note that e A %V ■ (a <8> a) is an analytic function, given by 

e A xPV- (o®o)(x) = (xPV- (a<g>a),</i(x- ■)), 

where gi the standard gaussian and the (•, •) refers to the duality product between com- 
pactly supported distributions and C°° functions. The properties of compactly supported 
distributions guarantee the existence of a compact K in M. d and C > , M G N such that 

( X PV-(a®a), 5l (x--)) <C J] supa Q 5l (^-y)<C" 5 i(x/2) 

M<M^ 

for large enough |x| . In particular, e A ^V ■ (a a) = C(|x| -5 ) as |x| — > oo. 

Let us now apply the asymptotic formula for convolution integrals (15) with g = g± and 
/ = (1 — x)A. We obtained this formula under the assumption / G E™ , with < i? < d. 
Here we have, instead, / G E™ C E™ but it is easily checked that such formula remains valid, 
in this case, also when d = 3, with the same proof, since / is locally integrable. Applying this 
formula in the case m = 2, and using J g± = 1 and J" y gi(y) dy = 0, we get, for |x| — > oo, 

e A (l - x)i(x) = 51 * f(x) = f(x) + 0(|x|- 5 ) 
= ^(x) + 0(|x|- 5 ) 
= PV- (a®a){x) + 0(\x\- 5 ). 



Theorem 2 is now completely proved. 



□ 
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